Increased accuracy of ligand sensing by receptor internalization 
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Many types of cells can sense external ligand concentrations with cell-surface receptors at ex- 
tremely high accuracy. Interestingly, ligand-bound receptors are often internalized, a process also 
known as receptor-mediated endocytosis. While internalization is involved in a vast number of 
important functions for the life of a cell, it was recently also suggested to increase the accuracy 
of sensing ligand as the overcounting of the same ligand molecules is reduced. Here we show, by 
extending simple ligand-receptor models to out-of-equilibrium thermodynamics, that internalization 
increases the accuracy with which cells can measure ligand concentrations in the external environ- 
ment. Comparison with experimental rates of real receptors demonstrates that our model has indeed 
biological significance. 
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I. INTRODUCTION 

Biological cells can sense and respond to various chem- 
icals in their environment. However, the precision with 
which a cell can measure the concentration of a specific 
ligand is negatively affected by many sources of noise 
[sf. Most noticeable is external noise from the random 
arrival of ligand molecules at the cell-surface receptors 
by diffusion. Nonetheless several examples exist in which 
measurements are performed with surprisingly high ac- 
curacy. In bacterial chemotaxis, for instance, fast moving 
bacteria such as Escherichia coli can respond to changes 
in concentration as low as 3.2 nM Q. This value is re- 
markable, since cells have only about one second between 
"tumbles" to evaluate the ligand concentration Q. Fur- 
thermore, this concentration value corresponds to only 
about three ligand molecules in the volume of the cell, 
assumed to be one fenitoliter, suggesting single molecule 
detection. High accuracy is observed also in spatial sens- 
ing by single cell cukaryotic organisms. Best character- 
ized are the slime mold Dictyostelium discoideum, which 
is able to sense a concentration difference of 1 — 5% across 
the cell diameter Q , as well as Saccharomyces cerevisiae 
(budding yeast), able to orient growth in a gradient of a- 
pheromone mating factor down to estimated 1% receptor 
occupancy difference across the cell @. Spatial sensing 
is also efficiently performed by lymphocytes, neutrophils 
and other cells of the immune system [l3| , as well as by 
growing synaptic and tumor cells. 

Previously, the fundamental physical limits to the ac- 
curacy of sensing as set by ligand diffusion have been cal- 
culated [TT| - [l5| . Recent work based on simplifying mod- 
els indicates that, if a cell effectively acts as an absorber 
of ligand, the accuracy is significantly increased [l6j . 
Such an increase in accuracy can be explained with the 
fact that absorption prevents ligand molecules from un- 
binding the receptors. Hence, the same ligand molecule 
can only be counted once by a receptor, avoiding a source 
of measurement uncertainty. However, whether cells with 
realistic receptors can act as absorbers and increase the 




FIG. 1: Single receptor, immobile at position xq, binds and 
unbinds ligand with rates k+c and fc-, respectively. 



accuracy of sensing is unknown. 

Motivated by these observations, in this paper we 
analyze the role of receptor- mediated endocytosis, i.e. 
the internalization of cither bound or unbound recep- 
tors from the cell membrane into the cell interior, often 
observed in eukaryotic cells [l3, [3- Internalization of 
ligand-bound receptors effectively leads to the absorption 
of ligand molecules and is therefore expected to draw the 
cell nearer to the physical limit of the perfect absorber 
[l6l |. Using simple models for the ligand-receptor dy- 
namics, we find that the effect of receptor-mediated en- 
docytosis indeed increases the accuracy of sensing ligand 
concentration, if internalization of ligand competes with 
ligand unbinding. Comparison of our results to the avail- 
able literature of experimental rate constants shows that 
receptors often work in this limit, indicating biological 
relevance of our results. 

The paper is organized as follows: In Section II we 
review the results regarding the accuracy of sensing for 
a single, immobile receptor without internalization. In 
Section HI, we study the role of internalization using a 
model of ligand-receptor dynamics with internalization. 
While ligand-receptor binding and unbinding is described 
by equilibrium thermodynamics, as previously developed 
in Refs. [T^iHOIj internalization clearly introduces non- 
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equilibrium thermodynamics into the problem. We con- 
sider the limit of fast diffusion, i. e. when the coupling to 
diffusion of ligand can be neglected, as well as the general 
case both near and far from equilibrium. This general 
case deals with the reduction of rcbinding of previously 
bound ligand molecules by receptor internalization, and 
hence represents the main result of the work. In Section 
IV we analyze the results obtained in Section III and 
connect to receptors from the biological literature. We 
conclude with final comments and discussion. An Ap- 
pendix is devoted to an alternative approach leading to 
the same results derived in the main text. 



II. REVIEW OF THE SINGLE RECEPTOR 

In this section we review previous results for a single, 
immobile receptor without internalization. Details of the 
method will be provided in Section III. As depicted in 
Fig. 1, such a receptor can bind and release ligand with 
rates fc+c and respectively. The kinetics for the oc- 
cupancy n{t) of the receptor are therefore given by 



dnjt) 
dt 



k+c[l ~ n(t)] - fc_n(t). 



(1) 



where the concentration of ligand c(z, t) = c is assumed 
uniform and constant. The steady-state solution for the 
receptor occupancy is given by 

(2) 
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with Ku = k-/k^ the ligand dissociation constant. The 
rates of binding and unbinding are related to the (nega- 
tive) free energy of binding through detailed balance 
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(3) 



with T the temperature in energy units. In the limit of 
very fast ligand diffusion, i.e. when a ligand molecule is 
immediately removed from the receptor after unbinding, 
the dynamics of the receptor is effectively decoupled from 
the diffusion of ligand molecules and hence, diffusion does 
not need to be included explicitly. 

Following Bialek and Setayeshgar [l^l , the accuracy of 
sensing is obtained by applying the Fluctuation Dissipa- 
tion Theorem (FDT) [1^, which relates the spectrum of 
the fluctuations in occupancy to the linear response to 
a perturbation in the receptor binding energy. Further- 
more at equilibrium the fluctuations in occupancy can be 
directly related to the uncertainty in ligand concentration 
using Eq. For a measurement performed on a time 
scale T much larger than the correlation time of the bind- 
ing and unbinding events, the fluctuations of the occu- 
pancy {{Sn)^) are obtained from the zero- frequency spec- 
trum divided by r. Using {{Sn)'^), one then obtains the 
uncertainty in measuring ligand concentration c [l2l . [20| 
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(4) 




FIG. 2: Receptor with internalization. In addition to ligand 
binding and unbinding described in Fig. 1, the receptor is 
internalized at rate ki, if occupied, and at rate fc^, if unoc- 
cupied. The internaUzed hgand is degraded. An unoccupied 
receptor is dehvered to membrane with rate kr. 



where the right-hand side is obtained for diffusion-limited 
binding (20j . i.e. when A:-|_c(I— n) — )■ 47rcD3S, with D3 the 
diffusion constant and s the dimension of the (spherical) 
receptor. Eq. @ shows that the accuracy of sensing is 
limited by the random binding and unbinding of ligand. 

In the case where diffusion of ligand is slow, ligand 
binding to the receptor is affected by diffusion The 
kinetics of the receptor occupancy and ligand concentra- 
tion are described by 



dn{t) _ 

dt 
dc{x, t) 
dt 



k-^-c{xo, t)[l — n(t)] — k^n{t) 



D3V^c{x,t) - d{x-XQ) 



dnjt) 



(5a) 
(5b) 



where xq indicates the position of the receptor and S{x — 
xq) is a Dirac delta function centered at the receptor 
location. The last term in the second equation describes 
a sink or source of ligand at xq , corresponding to ligand- 
receptor binding or unbinding, respectively. Analogous 
to fast diffusion, Eq. (5) has steady-state solutions c 
(independent of D3) and n given by Eq. (2). 

Following a similar procedure as in the previous case, 
the accuracy of sensing is given by [l^, [l^l 
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/s-l_c(l — n)T nsDjCT 
3 



2Tr SD3CT 



(6a) 
(6b) 



where the first term on the right-hand side of Eq. ((6aj) 
is the same as in Eq. Q, while the second term is the 
increase in uncertainty due to diffusion. This term ac- 
counts for the additional measurement uncertainty from 
rebinding of previously bound ligand to the rece ptor . For 
diffusion-limited binding, one obtains Eq. (|6b)) j2G | . 

Comparison of Eqs. (|4|) and (|6]) shows that removal 
of previously bound ligand by fast diffusion increases the 
accuracy of sensing, since the same ligand molecule is 
never measured more than once. 
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occupancy is given by 
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FIG. 3: Single receptor with internalization. Similar to Fig. 
[2l except that the delivery of receptor to the membrane is 
adjusted, so that Eq. ^ is fulfilled. This corresponds to an 
instantaneous replacement (round arrow) of the internalized 
receptor by a new receptor via the endocytosis machinery. 



III. EFFECT OF RECEPTOR 
INTERNALIZATION 

Here, we consider the case of receptor internalization. 
As depicted in Fig. 2, receptors at x = can bind and 
unbind ligand with given rates. Furthermore, a bound 
receptorStren can be internalized at rate fc^, while an 
unbound receptor can be internalized at rate kf. Hence, 
the kinetics of the fractions of occupied receptors n{t) 
and unoccupied receptors m{t) are given by 



dn{t) 



kj^^cm{t) — (fc_ + ki)n{t) 



(7a) 



^^^^ = -k+cm{t) - k^m{t) + k_n{t) + kr. (7b) 

Imposing a single receptor at x = a;o at any time via 

n{t) + m{t) = I, (8) 



Eq. (j7bp becomes redundant. As shown in Fig. [31 this 
condition implies that an internalized receptor is imme- 
diately replaced by a new, unoccupied receptor with rate 
kr{t) = k^m{t) + kin{t). Furthermore, rate fc^o of incom- 
ing ligand compensates for internalized ligand. 



A. Limit of fast diffusion 

We first consider the case of fast diffusion, i.e. when 
ligand unbound from the receptor is immediately re- 
moved. In this case the kinetics for the occupancy n(t) 
of the single receptor is described by 



dnjt) 

at 



= k+c[l - n{t)] - (fc_ + k,)n{t) 
= fc^c[l — n{t)] — K_n{t), 



(9) 



c + K^/k^ c + K 



M 



(10) 



where Km = is a Michael- Menten-type constant, 

and subscript i is used to indicate the steady-state value 
for the occupancy of the receptor in presence of internal- 
ization (c/. Eq. (2)). 

While Eq. (9) could be solved immediately by analogy 
to Eq. (1), we adopt here, for the fast diffusion case, 
the method of the effective temperature, which allows us 
to solve the general case in section HI B. Similar to the 
equilibrium case, at the non-equilibrium steady state the 
rates can formally be related to the binding free energy 



kj^c 



k-a + t) 



1 + h/k^ 



where we introduced the effective temperature 

T 

T, = 



ln(l+fci/fc„) 
In {k^c/k- ) 



1 



(11) 



(12) 



Hence, the effective temperature maps the non- 
equilibrium steady state to an effective equilibrium, 
allowing the generalization of the FDT to out-of- 
cquilibrium phenomena l2l[-[24l with applications in mod- 
cling biological processes [25|. Conceptually, an effective 
temperature larger than the environment temperature 
(Te > T) corresponds to a decrease in the receptor oc- 
cupancy, approximately reflecting internalization in the 
non-equilibrium steady state. 

In order to calculate the spectrum of the fluctuations in 
receptor occupancy, we follow Refs. [H,!!^] and consider 
small fluctuations around the stationary solution 

n{t) = fii + 6n{t). 

In order to apply the generalized FDT (gFDT) , we intro- 
duce fluctuations in the conjugate variable of the receptor 
occupancy, i.e. the free energy F, via fluctuations of the 
binding and unbinding rates 



5F 



Sk+ 

"fc7 



K- 



(13) 



where we approximate as a parameter. 
Linearization of Eq. (jH) leads to 



d[Sn{t)] 
dt 



-(fc+c + K_)(5n(t) 4- fc+c(l — rii) 



(14) 



where we used Eq. ([13)) to replace the fluctuations in 
the rate constants with fluctuations in the free energy, as 
well as steady state solution Ec^. (10). 

Fourier Transforming Eq. yields the susceptibility 



where k_ = k-+ki is the combined rate of unbinding and 
internalization. The steady-state solution for the receptor 
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describing the linear response of the receptor occupancy 
to a perturbation in the free energy. We now use the 
gFDT to calculate the spectrum S'„(w) ~ {\Sn{uj)\'^) of 
the fluctuations in n{t) 



5'„(w) 



2T, 2fc+c(l-n,) 
Im x(cj) = -- — z 77; o 



2{{5n) 



1 + (a;rc)2 ' 



(16) 



where the correlation time tc — ik+c + «-) ^ and the 
total variance 



27r 



fc+C + K- 



(17) 

have been introduced. Using Eq. ([TUl) . we calculate the 
uncertainty in ligand concentration from fluctuations in 
occupancy 



5c 



{c + Km)\__ Km 



-Sn = 



Km (l-ni)^"'" rij(l-ni) 

From Eq. (jlSp the normalized variance can be obtained 
{(Sc?) 1 



ni{l ~ Hi) ' 



(18) 
ned 

(19) 



corresponding to an instantaneous measurement. 

In the more realistic case, in which a measurement is 
performed during an averaging time r 3> tc, the error 
in the occupancy is linked to the low frequency spectrum 
via 



T k+CT 



(20) 



Using Eqs. ((T8| and <\2Q\ . the accuracy of sensing is given 
by 



{[5cf)r 



1 



fc+c(l — ni)T 27rZ?3Csr 



(21) 



Eq. ([21]) is identical to the result in Eq. (jlj without 
internalization, except that fii < n due to internalization. 
In fact, the removal of unbound ligand by fast diffusion at 
equilibrium is equivalent to removal of bound ligand by 
internalization at the non-equilibrium steady state. This 
equivalence can be readily seen from Eq. which is 
indistinguishable from simple unbinding with rate k_ = 
k^ + ki. Hence, the effective temperature applied here is 
in fact exact. 



B. Solution near equilibrium 

When considering ligand diffusion, the above proce- 
dure still applies with the exception that the concen- 
tration of ligand is allowed to vary due to binding and 
unbinding. The kinetics of the receptor occupancy and 
ligand concentration is described by 



dn{t) 
dt 
dc{x, t) 
dt 



k^c{xi^,t)[l — n{t)] — K_n(t) 

'dn{t) 



dt 



(22a) 

kin{t) 
(22b) 



where used as before. Furthermore, a 

source of ligand with rate fcoo is considered at location 
Xoo so as to compensate the loss of ligand molecules due 
to internalization. 

We note that the steady-state solution for the con- 
centration is not spatially uniform, but is depleted near 
the receptor due to internalization. This leads to the 
anomaly that we mathematically evaluate the rate of 
binding using the ligand concentration at Xq in Eq. (|22ap , 
while physically the diffusive flux and hence binding of 
ligand is determined by the ligand concentration c far 
away from the receptor (see section III C). This is reme- 
died by linearizing the ligand concentration around c in 
the following. Furthermore, Eq. ([8|) is again assumed 
valid, and therefore an additional equation describing the 
unoccupied receptor fraction m(t) with rates k^ and kr 
is redundant. 

Linearizing Eqs. (|22ap and p2b[) leads to 



d[5n{t)] 
dt 

d[5c{x, t)] 
di 



= fc+(l — ni)5c{xt),t) — {k+c + K-)Sn{t) 



+ (5fc+(t)c(l — Hi) — niSn- 
-- DsV'^dcix, t) - 6{x - xo) 

-ki5n{t) 



(23a) 



d[5n{t)] 
dt 



(23b) 



By applying the quasi-equilibrium picture with the effec- 
tive temperature Tg introduced in the previous subsec- 
tion, we use Eq. p3)) to introduce fluctuations in the 
free energy, and obtain 



J 
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^^^^^^^ = ^+(1 - ni)5c{x(3, t) ~ {k+c + K^)Sii{t) + K^Ui^ 



d[Sc{x,t)] 
dt 



d[6n{t)] 
dt 



kiSn{t) 



(24a) 
(24b) 



Fourier Transforming Eqs. p4ap and ()24b[) . wc obtain 

(25) 



D^q^ — iuj 

which can be inverse-Fourier Transformed in xq 

1 



(5c(xo,i^) = {ioj — ki)5n{uj) 



(26) 



Inserting Eq. ([26)) in the Fourier- Transformed Eq. (|24aP , 
we obtain 



Sn{uj) 1 



SF{uj) Tq fc+c + K_ + {ki — ia;)Ei(cj) — iw' 

(27) 

where Si(w) is given by 



Si(c^) 



d^q k+{l~n{) 



dq q'^[k+{l - n,)] 



(27r)3 Daq^ - iuj 2Tr^ D^q^ - iuj 

fc+(l - 7li)A fc+(l - n,;) 



2'KD:iS 



(28) 



Here A ~ tt/s is a cut-off due to the finite size s of the re- 
ceptor, introduced to regularize the integral in Eq. (|28p . 
As before, we apply the gFDT to derive the spectrum of 
the fluctuations Sn{uj) 

2Te 

Sn{uj) = — ^Im[x(cj)] 

2k+c{l - n,)[l + T.i{u})] 
~ [A:,Ei(u;) + k+c + K-]^ + oj^l + Si(w))2 ' 

(29) 



where we used Eq. (|27p for the susceptibility. In the 
realistic case, in which the measurement is time averaged 
over duration t much larger then the correlation time of 
the fluctuations, the relevant part of the spectrum is the 
zero frequency limit 



Sniuj ~ 0) 



2fc+c(l-fl,)[l + Ei(0)] 

[fc+C-f K_ +/CiEl(0)]2 



2{{Sr 



(1 + Ei(0))tc 

(l+fc,Ei(0)Tc)2 



~ 2(((5n)2)[l + aEi(0)]Tc, (30) 

where a = 1 — 2/s,;Tc7 < 1, and higher order terms in Ei(0) 
are neglected for sufficiently fast diffusion. As before, 
(((5n)2) = nj(l — Ui) and tc = (fc+c+ k^)^^. 



Using Eq. ()18|) , the normalized variance of the concen- 
tration is given by 



[ni{l - nj)]^ [ni{l-ni)] 



, (31) 



where r is the averaging time. Using Eq. (pO|) for the 
power spectrum, wc finally obtain for the accuracy of 
sensing with ligand internalization and diffusion 



a 



k+c{\ — rii)r ttDjcst 
l + 2a 



2'kD^cst 



(32a) 
(32b) 



The following conclusions can be drawn by compari- 
son with the result Eq. (jB]) without internalization: (i) 
Receptor internalization mainly reduces the second term 
in Eq.®, demonstrating for the first time that inter- 
nalization increases the accuracy of sensing by reducing 
the uncertainty from rebinding of previously bound lig- 
and. The first term is only reduced by replacing n by 
fii, with fii < n. (ii) In the limit fci — >■ 0, a — > 1 and 
the equilibrium result Eq. (jS]) without internalization is 
recovered, (in) As Eq. ([5^ becomes unphysical in the 
limit oi ki 00, i.e. does not approach Eq. (|2T|) without 
rebinding, our result can only be regarded an approxima- 
tion valid near equilibrium. 



C. Comparison with perfect absorber 

The perfect absorber is here defined as a receptor, 
which internalizes a ligand immediately once it is bound. 
Following Rcf. (T^ , the accuracy of sensing can be calcu- 
lated from the Poisson statistics of the number of binding 
events N in time r 



((<5c)2). {{SNy 



4:71 D3SCT 



(33) 



obtained from the diffusive flux of ligand to an absorb- 
ing sphere of radius s. Comparison with Eq. (|2ip for 
diffusion-limited binding without rebinding shows that 
the perfect absorber is yet more accurate by a factor 2. 
This is due to the fact that the fluctuations in occupancy 
in Eq. (|2ip stem from the random binding and unbind- 
ing/internalization events, while the uncertainty in Eq. 
P3p solely stems from the random binding events (see 
also factor 2 in Eq. (|A2p . as well as Ref. [2a] for further 
explanation) . 
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IV. IMPLICATIONS FOR BIOLOGY 

In the previous section we showed that internahzation 
increases the accuracy of sensing by reducing the mea- 
surement uncertainty from rebinding of previously bound 
ligand. In this section we review some receptors of known 
rate constants. We specifically would like to determine 
if the rate of internalization ki is fast enough, i.e. com- 
parable to the unbinding rate fc„, in order to effectively 
increase the accuracy of sensing. 

In Table I, we summarize experimental values for rate 
constants, including internalization, of various receptors. 
Most G-protein coupled receptors (GPCR) undergo in- 
ternalization [l^ . The Ste2 receptor in haploid yeast cells 
of a-mating type is involved in a pheromone sensing and 
signal transduction, leading to cell polarization, 'shmoo' 
formation and mating. The folate receptor (FR) in Dic- 
tyostelium, likely a GPCR [l^, H^, is used to sense and 
hunt bacteria. (The folate-binding protein in mammalian 
cells is a diagnostic marker for various cancers, and its 
internalization is exploited for drug delivery into cancer- 
ous human cell jHf.) However, the cARl receptor in Dic- 



tyostelium, used for sensing of cAMP under starvation, is 
not internalized [soj . The epidermal growth factor recep- 
tor (EGFR), a tyrosine kinase, is another important ex- 
ample of a receptor which is internalized [sij . This recep- 
tor is involved in cell growth, proliferation and differen- 
tiation (31I433I . Another class of internalized receptors is 
involved in uptake. Transferrin receptor (TfR) is used for 
iron uptake from extracellular space and plays therefore 
an important role in blood cells [31[. Transferrin binds 



to TfR, is internalized and releases its iron load through 
ion pump-induced pH reduction. The ligand-bound TfR 
is then recycled back to cell surface. Another example is 
the low density lipoprotein receptor (LDLR) [3l| . When 
bound to LDL-cholesterol via adaptin, LDLR is inter- 
nalized via clathrin-coated vesicles |3l| . Furthermore, 
the vitellogenin receptor (VtgR) is involved in oogenesis 
(egg- formation) [3l|. Once internalized, vitellogenin is 
turned into yolk proteins. Ligand-free receptors are re- 
cycled back to cell surface. Table I shows that in most 
cases ki is of the same order or larger than except 
for FR, where internalization is much slower than the 
unbinding of the ligand from the receptor. 



Receptor 


Function 


fc_ (min ^) 


Kd (nM ) 


ki (min '') 


fc^ (min"") 




chemotaxis: 










FR 


feeding 


0.096 [34] 


20.0 [34]t 


9.6x10"* [29]* 




Ste2 


mating 


0.06 [35] 


22.1 [35] 


0.24 [36] 


0.024 [36] 






0.0108 [32 


6.0 [37] 


0.156 [38] 


0.0156 [38] 


EGFR [31] 


development 


0.24 


2.47 


0.15 


0.02 




uptake: 










TfR [31J 


iron 


0.09 


29.8 


0.6 


0.6 


LDLR [31J 


cholesterol 


0.04 


14.3 


0.195 


0.195 


VtgR[3^ 


vitellogenin 


0.07 


1300 


0.108 


0.108 



TABLE I: Summary of experimental data for relevant receptor rates discussed in the main text. ' However, other values have 
been reported as well [4l|].* This rate is measured for folate-binding protein in cancerous mice cells, not in Dictyostelium. 



Figure 4 visualizes the contribution of internalization 
to the accuracy of sensing for the receptors from Table 
I. The faster the internalization, the larger the increase 
of the accuracy of sensing. However internalization can 
only reduce the second term in Eq. (j6ap from rebind- 
ing of previously bound ligand, not the first term from 
random binding and unbinding. To illustrate the rela- 
tive contribution of the two terms, we plot the square 
root of their ratio in Fig. 4 (filled bars). This shows 
that in most cases, in which internalization occurs, the 
noise ratio of the two terms is significant, ranging from 
few hundredths to order of unity. Hence, internalization 
can lead to a substantial increase in the accuracy of sens- 
ing. However, measured rate constants are substantially 
uncertain (see below) and diffusion coefficients of small 



ligand molecules range from 0.1/^m^/s in the syna ptic 
cleft between neurons I39l to 1 — 10//m^/s in blood [40| 
to 300^m^/s in water [l6|. In order to avoid this un- 
certainty in parameters, we also plot the upper limit of 
the noise ratio for diffusion-limited binding equal to \/2 
(dashed line in Fig. 4). Removal of the second noise 
term in this limit by internalization would increase the 
accuracy in Eq. (6b) by a factor 3. In Fig. 4, we also 
show the strength of internalization, defined by the ratio 
of internalization and unbinding rates (open bars) . 

How reliable are the measured values for the rate con- 
stants? Rate constants are generally obtained through 
radioactive labeling of ligand, with the isotope choice 
targeted to each specific case (the isotope ^^^/ giving 
the most accurate measurements) [43 |. In order to mea- 
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10 



- I -3 



- 0.1 



FIG. 4: Analysis of receptor rate data from Table I. (Filled 
bars) Noise ratio defined as the square root of the second (re- 
binding) term and the first (random binding and unbinding) 
in Eq. (6a). (Open bars) Strength of internalization, defined 
by ratio fci/fe-.The dashed line indicates the noise ratio {V2) 
for diffusion-limited binding (fc+c(l — n) — >■ 4ncDss) . Nu- 
merical values of the plot are provided in Appendix B. 



sure the unbinding and the binding rates, i.e. k_ and 
= k^/Ku, the receptors on the membrane must be 
separated by filtration or centrifugation from the soluble 
ligand. If the ligand-unbinding process is slow compared 
to separation, then the measurement of the amount of 
bound ligand through the radioactive label can be easily 
carried out; in the case of fast unbinding, measurements 
are less accurate. For the internalization rate, the ratio 
between the intensity at the surface and inside the cell is 
measured and from the slope of the time variation of this 
ratio, ki is determined |43|. This method, though, does 
not take into account recycling (diacytosis) of the recep- 
tor, which is a recurrent feature of the internalization 
process, also included in our model. These shortcom- 
ings, as well as the variability associated with different 
cell preparations, lead to a large error in the rate deter- 
mination 20%) [13 and variability between different 
measurements (20-90%). Other measurement methods 
employed in experiments include protease sensitivity as- 
says 44 [ and destination assays [45| . 

Many examples in fact suggest a direct relation be- 
tween ligand internalization and the accuracy of sens- 
ing, measured by the sensitivity of cell polarization or 
cell movement in shallow chemical gradients. These in- 
clude sensing of a- factor by budding yeast [13, IssI, fo- 
late by Dictyostelium [53] and PDGF by fibroblasts [46} . 
Other examples relate to embryonic development. In ze- 
brafish, primordial germ cells migrate toward chemokine 
SDF-la that binds and activates the receptor CXCR4b. 
It was recently shown that ligand-induced CXCR4b in- 
ternalization is required for precise arrival of germ cells 
at their target destination |47{ . During Drosophila oo- 
genesis, border cells perform directional migration [3^ . 
EGFR, together with two other receptor tyrosine kinases. 



is the main guidance receptor. Recent work in this sys- 
tem provided compelling evidence that guided cell move- 
ment also requires s pat ial control of signaling events by 
endocytic dynamics |48l . l49j . 



V. DISCUSSION AND CONCLUSIONS 

In this paper we analyzed the role of receptor internal- 
ization in the accuracy of sensing ligand concentration. 
By extending equilibrium single receptor models to non- 
equilibrium thermodynamics introduced by internaliza- 
tion, we derived expressions for the uncertainty in sens- 
ing ligand concentration. As expected, internalization 
of ligand-bound receptors makes the cell act similarly to 
an absorber and increases the accuracy of sensing. We 
then analyzed relevant experimental data, summarized 
in Table I, and concluded that in most cases, the con- 
tribution of receptor internalization to the increase in 
accuracy of sensing is non-negligible. However, a per- 
fect absorber is yet more accurate as its uncertainty only 
stems from random binding events, not from additional 
random unbinding and internalization events. Whether 
cells have developed mechanisms to approach the limit 
of the perfect absorber, e.g. by internalization, is not 
clear yet. However, since the accuracy of concentration 
sensing can always be improved by increasing the av- 
eraging time, one might expect that receptor internal- 
ization becomes increasingly important when time is of 
the essence. In addition to chemotaxis, embryonic de- 
velopment may be a biological system for which receptor 
internalization is important. Specifically, the Fgf8 mor- 
phogen, which regulates tissue differentiation and mor- 
phogenesis in zebrafish, forms exponential concentration 
gradients by diffusion and degradation, the latter being 
achieved precisely by receptor internalization [sof . How- 
ever, based on our results, internalization may also be 
used for the accurate readout of the gradient in the short 
amount of time dictated by cell division. 

Cells generally have many receptors to estimate ex- 
ternal concentrations of chemicals, leading to a spatial 
averaging and consequently further increase in the accu- 
racy of sensing. However, even the employment of many 
receptors cannot increase the accuracy of sensing beyond 
the physical limit of the perfect absorber [26j. In fact a 
cell only needs a relatively small number of receptors to 
achieve an accuracy comparable to the physical limit [3] . 
On the downside, if a cell uses many receptors, it needs 
to integrate this information in signaling pathways. If 
this process is fundamentally limited by noise as well, 
then this noise provides an upper limit on the overall 
estimation performance. 

Although we analyzed the role of receptor internaliza- 
tion in increasing the accuracy of sensing, it is important 
to stress that internalization fulfills several other known 
purposes in cells [sij . Among them are (i) redistribution 
of receptors to different locations on the cell membrane, 
(a) uptake of nutrients and chemicals, (in) signaling by 
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ligand in cell interior, and (iv) turning off persistent sig- 
nal as part of adaptation. All these aspects would need 
to be considered to fully characterize the working of a 
receptor. 

In order to derive the accuracy of sensing with internal- 
ization, we made a number of simplifying assumptions. 
In our model we neglected other possible sources of fluc- 
tuations such as fluctuations in receptor density, in order 
to relate our results to the single, immobile receptor. Fur- 
thermore, we introduced the effective temperature Te to 
generalize the FDT to non-equilibrium thermodynamics. 
While Te is well defined for well separated time scales, 
a potential time or frequency dependence of Tg [22| - [2^ 
was neglected here. However, our non-equilibrium result 
with internalization is consistent with the equilibrium re- 
sult in the fast diffusion limit. Removal of ligand by fast 
diffusion at equilibrium is equivalent to removal of bound 
ligand by internalization at the non-equilibrium steady 
state, providing confidence in our method. 

Internalization is not the only mechanism, by which a 
cell can act as an absorber and increase its accuracy of 
sensing. Other potential mechanisms include enzymatic 
degradation of ligand at the cell surface, e.g. degra- 
dation of cAMP by mPDE in Dictyostelium and of a- 
mating pheromone by Barl in budding yeast [13, l53j . 
Furthermore, at excitatory neural synapses, fast diffusion 
of AMPA receptors on the post-synaptic membrane sur- 
face has an important role in the sensing of neurotrans- 
mitter glutamate [s^ ]. Ligand-bound, desensitized re- 
ceptors diffuse away and are replaced by fresh receptors, 
leading to fast recovery and readiness for the next action 
potential and release of neurotransmitter. By the same 
mechanism, the accuracy of sensing may be increased, 
since ligand-bound receptors diffuse away and release lig- 
and far away from region of signaling, th us p reventing an 
overcounting of same ligand molecules. |55j . 
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Appendix A: Langevin approach 

As an alternative derivation, here we provide the solu- 
tion for the accuracy of sensing for internalization using a 
Langevin approach. We first consider fast diffusion. We 
start from Eq. © but add a noise term, ^„(t) 



dnjt) 



k+c[l - n{t)] - K_n{t) + ^„(i). 



(Al) 



where we assume 

(||„(w)p) = k+c{l - n,) + K.n, = 2fc+c(l - n,) (A2) 

due to Poisson statistics [IB, [131 • Linearizing and Fourier 
Transforming Eq. (jAip . assuming the rates fc+ and k_ 
constant, leads to 



Sh{uj) 



fc+C + K_ 



(A3) 



Hence, the power spectrum of the fluctuations in receptor 
occupancy is given by 



5„H = {\6n{uj)f) 



(fc+C -|- -I- iu' 

2fc+c(l — fii) 



r, (A4) 



where in the last step the property Eq. ()A2[) was used. 
Eq. (|A4[) is indeed equivalent to result Eq. (|16p in the 
main text. 

For the general solution, we start from the Fourier- 
Transformed Eq. ([24]) and Eq. ([Ml), i-e. 



and 

(iw - ki) 



dc{xQ,uj) 



fc+(l - m) 



Ei{Lo)Sh{uj)+U^), (A6) 



respectively, where 5n(w) and Cc(w) are additive noise 
terms, and Ei(a;) is given by Eq. pS]). Inserting Eq. 
in Eq. (|A5p and solving for Sh{uj) leads to 



dh{u!) 



(A7) 



fc+c + K_ + ki'Si{u}) — iuj{l + T,i{u})) ' 
from which the following expression for ensues 

.|2^^ fci(l-n.)^(|e;HP) + (||nHP) 

(A8) 

In the limit uj 0, using Eq. (|A2I) as in the previous 
case, we obtain 



Following [s^, we set 



Si(0) 



D3S k+[l-ni) 



in Ec^. (jA9p . and obtain for the power spectrum 
2((<5n)2)(H-I]i(0))Tc 



[1 + A:,I]i(0)tc]2 



(A9) 



(AlO) 



(All) 



Eq. (jAlip is identical to result Eq. (PH]) . obtained with 
the gFDT in main text. 
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Appendix B: Numerical values 

In this section we provide the numerical values used 
for plotting Fig. 4. Noise ratio and internalization 
strengths are respectively given by: 0.1/0.01 (FR), 0.08/4 
(Ste2), 0.47/0.625 (EGFR), 0.08/6.6 (TfR), 0.08/4.875 
(LDLR), 0.011/1.54 (VtgR). Specifically, to calculate the 



noise ratio, the first and the second term in Eq. (6a) 
are given by (in units of r): 2500/27.6 (FR), 4000/25.1 
(Ste2), 1000/224 (EGFR), 2667/18.6 (TfR), 6000/38.72 
(LDLR), 3428/0.426 (VtgR). We have used s = Inm, 
71 = 1/2, z. e. setting c = Kjj from Table I, and 
D3 ~ I /im^/s . 
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